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Abstract 

We derive cutting rules for loop integrals containing propagators with complex masses. Using a 
field-theoretical model of a heavy vector boson interacting with a light fermion, we demonstrate that 
the complex-mass scheme respects unitarity order by order in a perturbative expansion provided 
that the renormalized coupling constant remains real. 
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I. INTRODUCTION 



When describing resonances in perturbation theory, one needs to take their finite widths 
into account. A framework of handling this problem is provided by the complex-mass 
scheme (CMS) This approach has proven successful in various applications within 

the standard mo del 0-E3I. Recently, the CMS has also been applied in chiral effective field 



theory [13H15 



In the framework of quantum field theory for unstable particles the usage of the CMS 
leads to complex-valued renormalized parameters. While the problems of unitarity and 



causality in field theories containing unstable particles were resolved long ago [16| , the issue 



of perturbative unitarity of the S-matrix in the CMS is still open |4|. Using the CMS, one 
does not change the bare Lagrangian and, therefore, unitarity is not violated in the complete 
theory. On the other hand, perturbation theory is based on an order-by-order approximation 
to exact results. Therefore, it is not obvious that the approximate expressions for the S'- 
matrix also satisfy the unitarity condition. In the present work we examine perturbative 
unitarity in a model of a heavy abelian vector field interacting with a light fermion. Keeping 
the renormalized coupling constant (the expansion parameter of the perturbation theory) 
as a real quantity, we derive the cutting rules for one-loop integrals involving propagators 
with complex masses and show that unitarity is satisfied up to higher-order corrections in 
the coupling constant. In full agreement with Ref. 16[, the S'-matrix connecting only stable 
states satisfies the unitarity condition. 



II. IMAGINARY PARTS OF LOOP INTEGRALS 

In this section we propose a method for deriving the imaginary parts of (one-loop) inte- 
grals involving propagators with complex masses. In the next section we use these results 
to demonstrate the perturbative unitarity at one-loop order in a model of a heavy vector 
boson interacting with a light fermion. 

Before discussing the propagation of unstable particles, let us recall a few properties of 
the Feynman propagator of a stable particle. In momentum space the Feynman propagator 
of a scalar particle with mass m and four-momentum p is given by 1 
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In Eq. (jTJ), e is a positively defined quantity and the limit e — > + is assumed. Defining 
E p = a/p 2 + m 2 , the Feynman propagator has two simple poles at po = E p — it and po = 



i7i5(p 2 -m 2 ) . (1) 
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Besides the Feynman propagator, let us introduce the advanced and retarded propagators, 
respectively, 
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where a p = sign(p ) = — sign(— p ) = —o_ p . Note that A^(p) = A R (— p). In terms of these 
propagators, the Feynman propagator can be written as 



a f ( p ) = e( Po )A R (p) + e(-po)A A (p), 



(5) 



with the Heaviside step function 0(t) = 1 for t > and 6(t) = for t < 0. 

In the CMS, for an unstable scalar particle with mass M, width T, and four-momentum 
p, the Feynman propagator of Eq. ([1]) is replaced by 2 
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For a finite width T, the infinitesimal parameter e can be neglected. If one is also interested 
in the case of vanishing T, i.e. stable particles, one has to keep the infinitesimal parameter 
e. In the following, we drop e with the understanding that it is easily reintroduced by 
replacing M T — > M T + e. Let us define auxiliary functions which we denote as " advanced" 
and "retarded" propagators and corresponding to A': 
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2 We stress that in our notation the prime does not refer to the full, dressed propagator of the unstable 
particle. 
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where 
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Here and below it is understood that r<M and hence the small expansion parameter is 
r/M. The "advanced" ("retarded") propagator has two simple poles in the upper (lower) 
complex half plane and approaches the advanced (retarded) propagator of a stable particle 
as r 0. 

With the above definitions, let us consider a generic one-loop integral involving the prop- 
agation of both a stable and an unstable particle, 
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Using Eqs. ([!]) and ([6]), the imaginary part of the integral I\ is given by 
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The purpose of the subsequent manipulations is to bring the first term on the right-hand 
side of Eq. pip into a more convenient form. To that end, we make use of the following 
observation. As functions of the complex variable qo, both A^(q — p) and A' R (q) have 
simple poles in the lower half plane which is also true for their product. Closing the contour 
integration in the upper half plane including a vanishing contribution resulting from the 
semi circle at infinity, we find, using Cauchy's theorem, 
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By substituting Eqs. ([3]) and flHj) into Eq. ffl2|) and taking the imaginary parts of both sides, 
we obtain 
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Expanding x(q) — £ q + 0(T 2 ), we can replace Eq. ( TT3I) by 
d A k d 4 q 
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Subtracting Eq. (fill) from Eq. (fill and restoring e, we obtain 
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It is convenient to rewrite Eq. ( 1151) in a form analogous to the standard cutting formula for 
loop integrals with real masses [18|, ll9| . The left-hand side of Eq. (fT5l) is manifestly Lorentz 
invariant, and hence the right-hand side is Lorentz invariant order by order in T. Based on 
this fact we derive that 
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The symbol 0(T 2 ) on the right-hand side of Eq. ( 1T6|) indicates that the neglected terms 
contain at least one additional overall factor of V as compared to the first term. Note that 
the integrand of Eq. ( fl6l) is not obtained by expanding the integrand of Eq. ( 1151) . The 
equivalence of expressions ( 1T5"|) and ( TIB"]) can rather be seen by subtracting them from each 
other and considering, for e = 0, the limit r — > of the integrated result. Both Eq. ( I16p and 
Eq. (fT5|) turn into the expression for the stable particles in this limit. Hence the difference 
has to be suppressed by an additional factor of V. 

The generalization of the above procedure to any one-loop integral containing propagators 
with complex masses is straightforward. In analogy to Eq. (Fl6|) . cutting an unstable-particle 
line results in an overall factor of T, whereas cutting a stable-particle line generates a delta 
function. For integrals containing propagators of stable particles only, the usual cutting 
rules apply. A list of the imaginary parts of one-loop integrals needed for the calculations 
of the next sections is given in the appendix. 



III. THE MODEL 

In the next section we demonstrate that using complex renormalized masses in unstable- 
particle propagators does not violate perturbative unitarity. To that end let us consider a 
model describing the interaction of an unstable vector boson (B) with a stable fermion (-0), 

1 M 2 

C = -- iVi^f + ~Y B ^ B * + ^ W ~ m °) ^ + 9o 7^o (17) 

where F ^ u = d^B 0u — d u B ^. The subscript indicates bare parameters and fields. The 
masses are chosen such that the vector boson can decay into a fermion- ant if ermion pair. 
Since the vector boson couples to a conserved vector current, the above model is on- mass- 
shell renormalizable 1201, i.e. leads to finite physical quantities after renormalizing the masses 



and the coupling constant. 
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Our perturbative approach to the considered model is based on the path integral formal- 
ism. The integration over classical fields corresponding to particles (stable as well as unsta- 
ble) is performed in the standard way, i.e., the Gaussian part is treated non-perturbatively 
and the rest perturbatively. For stable particles, perturbation theory based on the path in- 
tegral formalism is equivalent to perturbation theory based on the operator formalism in the 
Dirac interaction representation. On the other hand, the functional integral also allows to 
incorporate the unstable degrees of freedom, while the application of the interaction picture 
using a "free" Hamiltonian for unstable states is conceptually problematic. 

We perform the renormalization in two steps: first we get rid off the ultraviolet divergences 
by applying dimensional regularization in combination with the MS scheme [ItJ • We refrain 
from showing the corresponding counter terms explicitly (including those leading to the wave 
function renormalization). Next we express the renormalized masses of the MS scheme in 
terms of physical quantities — the poles of the dressed propagators — and substitute them back 
into the Lagrangian. This amounts to performing the following substitutions in Eq. ( TT7T) . 

Bg B^, 

-> iff, 

itlq — > m + Sm, 

Ml M 2 - i MT + 5z = z + 5z , 

9o -> 9, (18) 



resulting in 



£ £main ~t~ £ct j 



£main = ~ F^F^ + | B ^ + $ (%$ - m) if, + Q $ 

5z 

C ct = —BpB"-8mif>if>. (19) 

The main Lagrangian £ ma j n generates the propagators and the vector-boson fermion interac- 
tion vertex. The counter-term Lagrangian C ct is treated perturbatively in a loop expansion, 
i.e., we write the counter terms as 5m = YlT=i^ jk ^' m k anc ^ ^ z = YlT=i^ k ^ z k an d include 
them order by order in perturbative calculations. 

The undressed propagators of the fermion and the vector boson take the following form, 
respectively, 

S F (p) = l — -, (20) 

i> — m + i e 

S'^P) = - ^WV£ (21) 
p 2 — z 

and, for later usage, we parameterize the self energy of the vector boson as 

i n^(p) = i [U(p 2 ) 9fiu + U p (p 2 ) PflPu ] . (22) 

IV. PERTURBATIVE UNITARITY OF THE S-MATRIX 

Below we demonstrate that the unitarity condition for the forward-scattering amplitude, 

Im [Tu] = (2tt) H\P n - Pi)TZ%i, (23) 
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is satisfied at one-loop order in perturbation theory. In Eq. (123]) . the T-matrix element 
between initial and final four-momentum eigenstates is written as (f\T\i) = (2n) A 6 A (Pf — 
Pi)7}i- The imaginary parts of one- loop integrals needed in this section are given in the 
appendix. 

We start with the tree-order amplitude of the process f(pi, (J\) + f(p2, o-z) f(P3, 03) + 
/(p4, C4) shown in Fig. 1 a), 



iTa 



where p = p% + P2 = P3 + P4 with p 2 = s and our fermion states are normalized as 
(f{p',a')\f{p,a)) = 2E(p)(2tt) 3 5(p' — p)b a i a (and analogously for antifermion states); the 
Dirac spinors are normalized as u(p, a)u(p, a) = 2m = —v(p,a)v(p,a). Because of current 
conservation, using Eq. (|2~T|) . we obtain from Eq. 



(25) 
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From now on we consider forward scattering, i.e., p\ = ps, P2 = £>4, 0\ = 03, and 02 = 04. 
Renaming —> V^, and using V* = V* 1 *, the imaginary part of the forward-scattering 
tree-order amplitude obtained from Eq. (|25|) reads 



Im [Ta] 
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s — z 


= Im 
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(26) 

We omit wave function renormalization constants for external fermion lines as they do not 
contribute to the obtained relations at the given accuracy. 

The one-loop diagrams contributing to the // — > ff amplitude are shown in Fig. [TJ 
Let us begin with discussing the imaginary parts of diagrams c) and i) relevant for the 
calculation up to and including order T 2 . Throughout we exploit the fact that g ~ 
(hence also ~ VT) • The results of these two diagrams read 



1% = iV^S'^{p)5z ig ^S' p MV h 



(27) 



where Il° /3 (p) and bz\ denote the self energy and the counter term at first order in h, 
respectively. Using Eq. (122]) . we obtain from Eq. ([2j 



Ti+c — T% + Tc 



^n l(g) + fel 



[s-z) 2 

f 

nary part of the one-loop order amplitude for s 7^ M 2 is obtained from Eq. 
Im [%+c 



(28) 



with n^s) denoting the one-loop contribution to the function n(s) of Eq. fl22|) . The imagi- 

as 



Im[(n 1 ( a ) + Sz^js-z*) 2 } 

\2( a _-y*\2 V V 



,S — Z)*[S — Z* 



im[n 1 (.) + ^ 1 ] (r3) _ (29) 

(s-z)(s-z*) M V ; v ; 
In Eq. (!29j) we took into account that g ~ y/T and hence ~ vT, Ili(s) ~ T, and 5z\ ~ Y. 
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k) 



i) 





m 



n) 



FIG. 1: Tree and one- loop contributions to // scattering. Solid and curved lines correspond to 
fermions and vector bosons, respectively. Crosses refer to counter-term contributions. 

Let us now consider the contribution of the intermediate state consisting of one fermion 
(/) and one antifermion (/) ( "square" of the tree-order amplitude shown in Fig. [2]) to the 
right-hand side of the unitarity condition of Eq. (123]) : 

- (2tt) - P)Tf fT)t T uf)i 



(2vr)5 [q 2 - m 2 ] Q{q ) (2vr)<5 [(p - q) 2 - m 2 ] Q(p - q ) 



2" ^ J (2tt) 4 

xu(pi,a 1 )j v v(p 2 ,a 2 ) S Pv (p) v(p - q,a)ipu(q,o) 

Xu(q, cr)-f a v(p - q, a) S' atl (p) v(p 2 , a 2 )^f, J ,u(p 1 , a x ) 
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In the second-to- last step we made use of the projection operators over the positive and 
negative energy states. The last step is obtained by applying the cutting rules for stable 
particles to the one-loop self energy of the vector boson, 



nf(p) = ig 2 



d A q 
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Using l + a q a p - q = 2[O(q )Q(p — q ) + Q(—q )O(q —p )], the fact that p > 2m in Eq. ([30]), 
and Eq. (122]) . results in Eq. ( 130]) . Let us compare Eq. (130]) with the sum of Eqs. (126]) and 
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FIG. 2: Squared s-channel tree-order // scattering amplitude. Solid and curved lines correspond 
to fermions and vector bosons, respectively, and the star indicates the complex conjugate of the 
amplitude represented by the corresponding diagram. 




FIG. 3: Squared t-channel, s-channel ^-channel interference, t-channel s-channel interference tree- 
order amplitudes. Solid and curved lines correspond to fermions and vector bosons, respectively, 
and the star indicates the complex conjugate of the amplitude represented by the corresponding 
diagram. 

( |29l) . Taking into account that Im[<5^i] = MT + 0(h 2 ) [which follows from Eq. ( |T8|) ]. we see 
that the unitarity condition is satisfied up to and including order T 2 . 

In the following, we qualitatively discuss the imaginary parts of the remaining diagrams 
of Fig. [T] using the formulae in the appendix. The sum of diagrams b) and d) may be 
expanded in powers of T. The first term in this expansion is real, whereas at second or- 
der the respective imaginary parts of the two diagrams cancel each other. Therefore, the 
corresponding imaginary part of these two diagrams is of 0(T 3 ). Similarly, expanding the 
amplitude corresponding to diagram j), the first term in this expansion is of order g 4 and is 
real and the imaginary part starts contributing at 0(T 3 ). The imaginary parts of diagrams 
f), g), h), k), and 1) are all of higher order. This is the case because all these diagrams 
are proportional to g 4 and applying vertical cuts (relevant to forward scattering) in each 
case means cutting at least one unstable-particle line producing an additional factor of T. 
Finally, cutting two stable fermion lines in diagrams e), m), and n) generates the imaginary 
parts corresponding, up to higher-order terms, to the "square" of the tree-order diagrams 
shown in Fig. |3j On the other hand, all remaining possible cuts involve at least one unstable 
line which leads at least to 0(T 3 ). 

According to Ref . [16( , only stable asymptotic states contribute to the unitarity condition 
of the S'-matrix. In order to demonstrate that our cutting rules agree with this result, we 
consider the one-loop contribution to the forward-scattering amplitude of /// — > fff shown 
in Fig. H]a). Using Eq. ( ITS]) , it is easily seen that the imaginary part obtained by cutting 
the two lines of the loop in diagram a) is proportional to T. In Fig. H]b), it is schematically 
represented as the "square" of the tree-order diagram (modulo higher-order corrections). As 
the width T is generated by diagrams representing the decay of the vector boson into stable 
particles, it is clear that the imaginary part of the diagram in Fig. H] a) corresponds to the 
"square" of diagrams with stable particles only in external legs [see diagram c) in Fig. H]. 
Since the width T is calculated for an "on-mass- shell" vector boson, diagram c) contains 
also contributions corresponding to loop diagrams of higher order. Note that in the limit 
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x 




a) 



b) 




c) 



FIG. 4: /// scattering, a) one- loop contribution, b) "square" of the tree-order diagram, c) dia- 
grams with stable particles only in external legs. Solid and curved lines correspond to fermions 
and vector bosons, respectively, and the star indicates the complex conjugate of the amplitude 
represented by the corresponding diagram. 

of vanishing T (M < 2m) one is not allowed to drop e in Eq. ([TrJj) . Therefore, the result of 
diagram H]b) does not vanish. In fact, for T = the limit e — > leads to a delta function 
corresponding to the vector line, and we obtain the standard cutting rule for stable particles. 

V. CONCLUSIONS 

In this work we developed a procedure for deriving the imaginary parts of (one-loop) 
integrals involving propagators with complex masses. With the aid of this method, we 
demonstrated perturbative unitarity of the scattering amplitude within the complex-mass 
scheme at the one-loop level. This result was obtained under the assumption that the 
expansion parameter of perturbation theory (the renormalized coupling constant) remains 
real. Our results are in full agreement with the findings of Ref. [l6| that unstable states do 
not appear as asymptotic states and are therefore excluded from the unitarity condition. A 
generalization of cutting rules for unstable particles to higher orders of the loop expansion 
is straightforward. However, because of the non-trivial dependence of imaginary parts on T, 
the analysis of perturbative unitarity in higher orders will become more involved. 
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Appendix A: Imaginary parts of one-loop integrals 

In order to compactify the notation let us introduce the following abbreviations, 

d(k 2 ) = (k 2 -m 2 ) 2 + e 2 , 
D(k 2 ) = (k 2 — M 2 ) 2 + M 2 T 2 . 

We consider the integral 

h = l j (S 4 (2^ (27r)4 ^ (fc + l ~ Pl) (2?r)4(54(fc + q - P2) A ' {k) Af(9) Af{1) 

A'(k)A F (p 2 — k)A F (pi — k) . (Al) 
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To rewrite the expression of Eq. ( 1A2I) in a more convenient form we consider the following 
integral, 
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take its imaginary part and obtain 
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By subtracting Eq. (IA4I) from Eq. (IA2I) . we obtain 
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The first term in Eq. (1A5I) corresponds to cutting both stable-particle lines [18l . Il9j . The 
other two terms correspond to cutting one of the two stable-particle lines together with the 
unstable-particle line. These two terms are proportional to T. 
Next, let us consider an integral 

d*k d*q dH d*p {2 ^ {k+px _ pz _ q){2 ^ {k+px _ l) 



(2tt) 4 (2tt) 4 (2tt) 4 (2tt) 4 
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= i J -0^ A'(k)A'( Pl -p 3 + k)A F ( Pl + k)A F (p 2 - k) . (A6) 
Its imaginary part reads 
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MT MT „, j2 2 N r/ 2 2n1 /a-n 

tt5(Z 2 - m 2 )7r5(p 2 - m 2 U . (A7) 



To rewrite the expression of Eq. (1A7I) in a more convenient form we consider the following 
integral, 

f r\ 4 k d 4 n d 4 l d 4 n 

= i t 1 1 l L 1 ^ (^) 4 S\k + Pl - P3 -q) {2*)W(k + Pl - Q 



(2tt) 4 (2tt) 4 (2tt) 4 (2vr) 4 
x(27r) 4 5 4 (fc -p 2 +p)A' A (k)A' A (q)A A (l)A R (p) 

= * f 7^4 A^)A^( Pl - p 3 + fc)A^(pi + fc)A fl (p2 - k) , (A8) 
take its imaginary part and obtain 

/d 4 k d 4 n d 4 l d 4 n 
(2tt) 4 (2tt) 4 (2tt) 4 (2tt) 4 v ; v F1 Fi yM ; v F1 ; 

k 2 - M 2 q 2 - M 2 I 2 -m 2 p 2 - m 2 



x(27r) i 5\k-p 2 + p) 



D(k 2 ) D(q 2 ) d(l 2 ) d(p 2 ) 



k 2 -M 2 q 2 -M 2 2 2 2 2 2 , Mr|Mr|;2_ m 2 2_ m 2 

■ n 0(1 — m ) d(p — m ) ai<j p 



D(k 2 ) D(q 2 ) v 7 v ^ ' 1 p D(k 2 ) D(q 2 ) d(l 2 ) d(p 2 ) 

MT I 2 - _ M 2 n 2 - m 2 MT¥ a 2 - M 2 I 2 - m 2 

1V1 £ k Q ivi p m 2 2\ , s k Q. ivi i m 2 2 



D(k 2 ) D(q 2 ) d(p 2 ) v y 1 D(£; 2 ) D(g 2 ) d(Z 2 ) 
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■nb(l -m)ai+ —. . —. . . . 7rd(p -m)cr p 



D(k 2 ) D(q 2 ) d(p 2 ) K 11 D(k 2 ) D(q 2 ) d(l 2 ) 

MT ^ MY f- 1 

fg 7T <5(/ 2 - m 2 ) 7T 5(p 2 - m 2 ) <t, a p + C(r 2 ) \ . (A9) 



D(fc 2 ) L>(g 2 ) 
Subtracting Eq. (1A9h from Eq. (I A 71) . we obtain 



x (2tt) V(A; - p 2 + p) 1-^77^ 4^ ^ " ^) W - ^X 1 + * °p) 



k D(g 2 ) 
Mr (i-lf) Mr Z 2 -m 2 p 2 -m 2 



D(/c 2 ) D(g 2 ) rf(/ 2 ) d(p 2 ) 

it ai ) q 2 — M 2 p 2 — m 2 



MT n 2_ M 2 n 2_ rr ,2 

n5(l 2 -m 2 ) 



D(k 2 ) D(q 2 ) d(p 2 ' 

MT ( 1 + t>) g 2 -M 2 Z 2 -m 2 



D(k 2 ) D(q 2 ) d(l 



2\ 



it 5(p — m ) 



£>(fc 2 ) D(g 2 ) ci(p 2 



7r o(/ -m J 



; " 2 7T(5(p 2 -m 2 )+ (9(r 2 )[>. (A10) 



£>(fc 2 ) D(q 2 ) d(V 

Note that we have kept the second line in the brackets for completeness although it is of 
OCT 2 ). In analogy to the step from Eq. (fT5l) to Eq. (flEI) . one can replace expressions like f 2 - 
by the corresponding functions a q in Eqs. (1A5j) and (lAlOf) . Finally, the case T = is obtained 
by first replacing M T — > M F + e in all integrals above. Setting T = and then taking the 
limit e — > 0+, we exactly reproduce the standard cutting formulas for loop integrals with 
real masses Fill, 19| . 



[1] R. G. Stuart, in Z° Physics, ed. J. Tran Thanh Van (Editions Frontiers, Gif-sur-Yvette, 1990), 
p. 41. 

[2] A. Denner, S. Dittmaier, M. Roth, and D. Wackeroth, Nucl. Phys. B560, 33 (1999). 
[3] A. Denner, S. Dittmaier, M. Roth, and L. H. Wieders, Nucl. Phys. B724, 247 (2005). 
[4] A. Denner and S. Dittmaier, Nucl. Phys. Proc. Suppl. 160, 22 (2006). 

[5] A. Denner, S. Dittmaier, M. Roth, and L. H. Wieders, Nucl. Phys. Proc. Suppl. 157, 68 
(2006). 

[6] A. Bredenstein, A. Denner, S. Dittmaier, and M. M. Weber, Nucl. Phys. Proc. Suppl. 160, 
131 (2006). 

[7] A. Bredenstein, A. Denner, S. Dittmaier, and M. M. Weber, JHEP 0702, 080 (2007). 
[8] S. Actis and G. Passarino, Nucl. Phys. B777, 100 (2007). 



13 



[9] S. Actis, G. Passarino, C. Sturm, and S. Uccirati, Phys. Lett. B 669, 62 (2008). 
[10] S. Actis, G. Passarino, C. Sturm, and S. Uccirati, Phys. Lett. B 670, 12 (2008). 
[11] A. Denner, S. Dittmaier, T. Kasprzik, and A. Muck, JHEP 0908, 075 (2009). 
[12] A. Denner, S. Dittmaier, T. Gehrmann, and C. Kurz, Nucl. Phys. B836, 37 (2010). 
[13] D. Djukanovic, J. Gegelia, A. Keller, and S. Scherer, Phys. Lett. B 680, 235 (2009). 
[14] D. Djukanovic, J. Gegelia, and S. Scherer, Phys. Lett. B 690, 123 (2010). 
[15] T. Bauer, J. Gegelia, and S. Scherer, Phys. Lett. B 715, 234 (2012). 
[16] M. J. G. Veltman, Physica 29, 186 (1963). 

[17] J. C. Collins, Renormalization (Cambridge University Press, Cambridge, England, 1984). 
[18] R. E. Cutkosky, J. Math. Phys. 1, 429 (1960). 

[19] M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory (Addison- Wesley, 

Reading, USA, 1995). 
[20] D. G. Boulware, Annals Phys. 56, 140 (1970). 



14 



